s18dcf

s18dcf © Numerical Algorithms Group, 2002.

Purpose

S18DCF Modified Bessel functions K (z), real a >= 0, nu + a complex z, nu = 0,1,2,...

Synopsis

[cy,nz,ifail] = s18dcf(fnu,z,n<,scale,ifail>)

Description

 
 This subroutine evaluates a sequence of values for the modified 
 Bessel function K    (z), where z is complex, -(pi) < arg z <= 
                  (nu)                                         
 (pi), and (nu) is the real, non-negative order. The N-member 
 sequence is generated for orders (nu), (nu)+1,...,(nu)+N-1. 
                                                   z
 Optionally, the sequence is scaled by the factor e .
 
 Note: although the routine may not be called with (nu) less than 
 zero, for negative orders the formula K     (z)=K    (z) may be 
                                        -(nu)     (nu)          
 used.
 
 For very large |z| or ((nu)+N-1), argument reduction will cause 
 total loss of accuracy, and so no computation is performed. For 
 slightly smaller |z| or ((nu)+N-1), the computation is performed 
 but results are accurate to less than half of machine precision. 
 If |z| is very small, near the machine underflow threshold, or 
 ((nu)+N-1) is too large, there is a risk of overflow and so no 
 computation is performed. In all the above cases, a warning is 
 given by the routine.
 

Parameters

s18dcf

Required Input Arguments:

fnu                                   real
z                                     complex
n                                     integer

Optional Input Arguments:                       <Default>

scale (1)                             string   'u'
ifail                                 integer  -1

Output Arguments:

cy (n)                                complex
nz                                    integer
ifail                                 integer